A Kolmogorov-Zakharov Spectrum in AdS Gravitational Collapse 
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We study black hole formation during the gravitational collapse of a massless scalar field in 
asymptotically AdSo spacetimes for D = 4, 5. We conclude that spherically symmetric gravitational 
collapse in asymptotically AdS spaces is turbulent and characterized by a Kolmogorov-Zakharov 
spectrum. Namely, we find that after an initial period of weakly nonlinear evolution, there is a 
regime where the power spectrum of the Ricci scalar evolves as uj~ s with the frequency, w, and 
s ~ 1.61 which is tantalizingly close to Kolmogorov's 5/3. 



Introduction. 

Due to the AdS/CFT correspondence the question of 
instability of Anti-de-Sitter (AdS) spacetimes sits in the 
intersection of mathematical and numerical relativity, 
string theory, field theory and condensed matter physics. 
The stability of spacetimes is central to the study of gen- 
eral relativity as spacetime is not just a stage where inter- 
actions take place, rather, it is a dynamical participant. 
AdS is the maximally symmetric solution to the Einstein 
equations with a negative cosmological constant and one 
would expect it to be stable. It has recently been ar- 
gued that AdS is non-linearly unstable |l| to black hole 
formation. Namely, an arbitrarily small generic pertur- 
bation leads to the formation of a black hole. This is 
completely different from the situation in the other two 
maximally symmetric solutions to the vacuum Einstein 
equation - Minkowski (zero cosmological constant) and 
de Sitter (positive cosmological constant) - which are 
known to be stable. The feature responsible for this in- 
stability is the fact that AdS has a timelike boundary at 
spatial and null infinity. The mechanism proposed in [1] , 
called weakly turbulent, as responsible for this instability 
was resonant mode mixing that gives rise to diffusion of 
energy from low to high frequencies. 

The study of non-linear instability of AdS initiated in 
[H has been followed up in Q for higher dimensions and 
in Q for a different type of configuration in pure AdS 
but with angular dependence. Some discussions in the 
context of the relation to thermalization in the dual field 
theory have appeared in 0, 0] ■ 

In this paper we study in more details the turbulent 
mechanism conjectured in [l[. We implement a sugges- 
tion put forward in to apply the arsenal of meth- 
ods of dynamical systems to the gravitational collapse 
in asymptotically AdS spacetimes. One of the most cel- 
ebrated results in the theory of turbulent flows is Kol- 
mogorov's theory of 1941 [7], a master piece of dimen- 
sional analysis (for modern accounts see the monographs 
[H and Q). In Q we argued that the power spectrum 
of the Ricci scalar at the origin shows period-doubling, 
a sign of turbulent evolution; here we establish that, 



in the appropriate regime, the power spectrum has a 
Kolmogorov-Zakharov scaling with the frequency, ui, as 
cj~ s , with s ~ 1.61 which is tantalizingly close to Kol- 
mogorov's 5/3. 

Gravitational collapse in asymptotically D- dimensional 
AdS spacetimes 

We consider the dynamics of a massless scalar field, tp, 
in D = d + 1 dimensions minimally coupled to gravity 
with a negative cosmological constant A: 



S= d d 



xy/-g 



1 



16ttG 



(R- 



(1) 



where G is Newton's constant. We focused on spherically 
symmetric configurations described by the following met- 
ricQ, 
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where £ 2 — —d(d — 1)/2A, dQ 2 is the metric on the unit 
(d — l)-sphere. The functions A, S and the scalar field, ip, de- 
pend on (t, a;). Notice that the spatial domain is contained in 
the interval < x < tv/2. The AdS spacetime, which is the 
maximally symmetric solution to the vacuum Einstein equa- 
tions with a negative cosmological constant, A, corresponds 
to A = 1, 5 = and ip — 0. In this context it is the analog of 
the Minkowski spacetime for zero cosmological constant and 
dc Sitter spacetime for positive cosmological constant. 

By introducing the auxiliary variables $ = tp' and II = 
A~ 1 e s ip, where the overdots and primes denote derivatives 
with respect to t, x, respectively, the field equations read: 
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The first two equations are constraint equations, the third 
equation is a consequence of the definition of the auxil- 
iary variables, and the last is the Klein-Gordon equation 
5 MI/ V nidvtp) = 0. Hereafter, we assume units where 8nG = 
d-1. 

There is a natural mass function, m(x,t), in AdS space- 
times given by 



2m r 2 

r d-2 ' ~~p, 



(7) 



where the standard spherical coordinate r is related to x as 
r = lta,n(x/ 1), and one can show that, 



m(x,t) = (1 - A)- 



2ooB"(f)' ■ (8) 
This expression gives the total mass-energy inside a radius x 
at the instant t. The ADM mass of the system is obtained 
by evaluating the mass function asymptotically, or Madm = 
lim, c _ > , r «/2 m(x, t). As a typical feature of various spacetimes, 
including AdS, the ADM mass is a conserved quantity. 

Toward a dynamical system via spectral methods 

We integrate numerically the field equations (2)- (5) using 
the Galerkin- Collocation method [lCj. The central notion of 
any spectral method is to approximate a system of partial 
differential equations by a finite set of ordinary differential 
equations associated to unknown coefficients or to the grid 
values of certain quantities. This dynamical system approach 
proves very useful in describing the complex dynamics of the 
scalar field collapse in D-dimensional AdS spacetimes. 

The starting point is to provide an approximation for the 
metric functions in the form of series of appropriate basis 
functions: 



n a = ]T a k (t)Mv), *« = Y, 6 *(*)x*(i/)> (9) 

fc=0 fc=0 

M M 

A a = i + J2^m ( k A) (y), 5 a = ]T M*)V>i%). (io) 

fc=0 fc=0 

In the above expressions M, N are the truncation orders; 
a,k{t), bk(t), Ck(t),5k(i) are the unknown modes or coefficients 
that constitute the spectral representation of the metric func- 
tions. The basis functions i'k{y),Xk(y),i'k A) '{v)^k\y) de " 
pend on the spatial coordinate y = ix/nl — 1 which varies 
in the interval — 1 < y < 1. These basis functions satisfy 
the condition of regularity of the field equations at the origin 
(x — or y = —1) and at the spatial infinity (x = for/2 or 
y = 1). Near the origin it follows that [1[: 



n(i, x) = Ho(t) + 0(Sy 2 ), $(t, x) = 0(Sy), 



(11) 



A(t,x) = l + U (t) + O(Sy 2 ), S(t,x) = 0(Sy 2 ), (12) 

where 5y = y + 1. The conditions near the spatial infinity 
depend on d, 



n(i, x) = H d (t)8y d + 0(Sy d+2 ), S(t, x) = 0{8 y d ^), (13) 

A(t, x) = l + A d Sy d + 0(Sy d+2 ), 5(t, x) = S 00 + 0(Sy 2d ), 

(14) 



where 5y = y — 1. 

The substitution of the approximations given by Eqs. ([9]) 
and (| 10p into the field equations ©-© results in the residual 
equations. We impose that the residual equations vanish at 
certain points known as the collocation or grid points. As 
a consequence, the field equations are transformed into sets 
of ordinary differential equations for the modes (a,k(t),bj(t)), 
along with algebraic equations Ck = Ck(ai, bj, Hi, <6 m ) and 
Sk = 5k(ai,bj,Ili,&m), where II;, $ m are the grid values of 
the functions II, 3>, respectively. 

To evolve the spacetime we consider the initial data 0]: 
$(0, x) = 0, 11(0,3;) = eo exp (— tan 2 x/a 2 ), with a and eo 
as free parameters. In the spectral space this initial data is 
expressed as bj(0) = and 11(0, Xj) = J2k=o a k(0)i>k{xj), 
where Xj denotes the collocation points. From this last ex- 
pression, the initial modes a,k(Q) are obtained. With these ini- 
tial conditions, the dynamical system is numerically evolved 
and the relevant functions reconstructed afterwards. We will 
show elsewhere the full details of this procedure and the 
numerical tests that guarantee the accuracy of the integra- 
tion. Here, as an example of numerical accuracy, we show in 
Fig. ((TJ) the behavior of the relative error in the ADM mass 

\M ADM {t)-M AD M(t=0)\ where M ig defined be l ow equa _ 

tion ([8}. The truncation orders are from top to bottom: 10 
(gray), 20 (cyan), 30 (red) , 40 (blue) and 50 (gold). Note that 
the periodic structure in time (with period approximately it) 
is the result of the timelike boundary at spatial infinity typi- 
cal of AdS and has been pointed out in numerical simulations 
in [2,1,0. 
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FIG. 1. Relative error of the ADM mass as a function of time. 



The weakly nonlinear regime 

By decomposing the scalar field in eigenmodes of the AdS 
Laplacian, the analysis of T] gained important insight into 
the interactions among these modes during the gravitational 
evolution. It was shown analytically that for generic initial 
data there is no way to remove all the resonances; some secu- 
lar terms persist at higher order in perturbation theory, even 
after accounting for modulation of phase. Inspired by this 
analytical result, [l|] also showed numerically that the time 
scale after which the weakly nonlinear approximation breaks 
down is of the order tNL ~ 1/eo, where eo is the amplitude of 
the initial profile. 

Let us discuss, motivated by the results of [J], a conve- 
nient way of viewing the dynamics of a scalar field in AdS 
spacetimes 0]. In the weak field regime, j,l-4< 1, the 
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Einstein-Klein- Gordon equation in spectral space can be writ- 
ten as, 



Ofe(t)- 



t&n 2 (x(y)) 



[tan 2 (a;(y))]'afe(t)+J r (5, 1— A, a k , d k , d k , ..) = 0. 

(15) 



The first two terms represent the evolution of the linearized 
modes (see Wald [ll]]) of the scalar field that oscillate in well- 
defined frequencies w| = (d + 2j) 2 , j = 0, 1, 2, .... The term 
denoted by T contains the interaction with the gravitational 
sector. Therefore, we imagine our system as a set of non- 
linear ly interacting oscillators, and consequently the periodic 
motion can be understood as taken place on a torus in the 
abstract phase space spanned by the modes a k (t). Further 
numerical experiments indicated that if the term T is ini- 
tially small, its influence on the long term dynamics is to alter 
the natural frequencies of oscillation producing quasi-periodic 
motion. Then, the motion becomes unstable and the torus is 
eventually destroyed forming a strange attractor, which is a 
signature of chaos in the dynamical system. This means that 
the evolution of the scalar field in AdS is turbulent. This 
route to turbulence resembles the Ruelle-Takens scenario pj| 
where the initial amplitude eo plays the role of the control 
parameter. We presented evidence for such picture in [6J. 
Nonlinear regime: A Kolmogorov-Zakharov spectrum 

For a diffeomorphism invariant theory physical information 
is contained in coordinate-independent quantities; this is one 
of the most challenging aspects of applying the standard the- 
ory of dynamical systems to gravitational problems. We fol- 
low up on one of the suggestions presented in [g : to study the 
power spectrum of the scalar curvature. For the backgrounds 
we are considering the scalar curvature can be written as: 



R = l/l 2 [Acos{xy{-lf 



d(d-l)]. 



(16) 



We have found similar results using the power spectrum of 
the mass defined in equation [8] Mindful that results about 
the power spectrum might be affected by coordinate redefi- 
nitions but hopeful that the general theorems proving that 
relativistic chaos is coordinate invariant [l3|] extend to our 
context, we proceed to a detailed study of properties of the 
power spectrum of the Ricci scalar. 

Given the notion that for a given amplitude eo the weakly 
nonlinear regime lasts for a time interval proportional to 1/eo, 
by suitably choosing the amplitudes we can reliably separate 
the regimes of the scalar field evolution during collapse. In 
Fig. ((2| we present the power spectrum for two systems with 
amplitudes 10 -3 (red) and 2.5 (black) respectively. The small 
amplitude evolution is still in the weakly non-linear regime 
where the system is well-described by a set of weakly cou- 
pled oscillators yielding well-defined peaks at frequencies cor- 
responding to those of the eigenvalues of the Laplace operator 
in AdS. This is precisely what we see in the power spectrum 
of this signal. Clearly, for higher frequencies we already wit- 
ness a more complicated picture signalling that the oscillators 
are not free and, as discussed in [l(, and interact by transfer- 
ring energy to higher frequency modes. In the same figure 
we have plotted the power spectrum for a higher amplitude 
simulation in which the nonlinear regimes sets in earlier. The 
spectrum of the of the nonlinear evolution seems to follow a 
precise scaling and we have also plotted the -5/3 slope line 
that corresponds to the Kolmogorov spectrum. 




FIG. 2. Power spectrum for different amplitudes in gravi- 
tational collapse in AdS^. The linear regime (eo = 10~ 3 ) 
with well-defined frequencies (black line); a nonlinear regime 
(eo = 2.5) that approximately follows a power law and a 
straight line that is the Kolmogorov 5/3 scaling (red). 



It is well known that Kolmogorov scaling is a property of 
a steady-state system, that is, of systems subjected to an 
injection of energy at some scale and allowing a sink or energy 
at some other scale UQ. Therefore, the only way for us to 
view the appearance of this type of spectrum during collapse 
is as a regime in the full evolution of the black hole formation. 
Since the time scales of the different regimes is governed by 
the initial amplitude eo we have de facto created a steady-state 
situation by suitably choosing the values of the amplitude. 

We are, therefore, not describing a steady-state situation 
but a phase in the evolution of scalar gravitational collapse 
in AdS. We know that the ultimate fate in the evolution of 
the system is the formation of a black hole. However, we un- 
cover Kolmogorov scaling by judiciously choosing the ampli- 
tude and the time interval from our numerical simulations of 
gravitational collapse. We have, at least two important scales 
in the collapse problem: the nonlinear scale t^L ~ 1/eo an d 
a new scale that we call the Kolmogorov scale, tx\ evolution 
beyond tn involves the violent collapse into a black hole. 

In Fig. ((3| we present the power spectrum for the Ricci 
scalar of a simulation with amplitude eo = 3.78. The spec- 
trum for high frequencies is difficult to interpret. However, re- 
stricting to relatively low frequencies indicates a Kolmogorov- 
Zakharov spectrum. 

Turbulence is very difference depending on the dimension- 
ality of the space. Below we provide results for collapse in 
asymptotically AdS± spacetimes. However, in the spherically 
symmetric approximation to gravitational collapse all func- 
tions depend on only two coordinates - time and radius. As 
can be seen from equation different powers of the radius 
dominate the evolution near the center. We verified explicitly 
a claim advanced in 0] that the same turbulent mechanism is 
at works in spherically symmetric collapse in AdS 4 and AdSs . 
Below we present the power spectrum for gravitational col- 
lapse in AdSi and verify the existence of a power scaling in 
the appropriate regime see Fig [4] 
Conclusions 

We have provided a more complete picture of the gravita- 
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FIG. 3. The spectrum for black hole simulations with ampli- 
tude eo = 3.76. 
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FIG. 4. Power spectrum for different amplitudes in gravi- 
tational collapse in AdSi. The linear regime (eo = 10 -2 ) 
with well-defined frequencies (black line); a nonlinear regime 
(eo = 2.5) that approximately follows a power law and a 
straight line that is the Kolmogorov 5/3 scaling (red). 



tional collapse in AdS. In particular, we have gone beyond 
previous descriptions which focused on the weakly nonlinear 
regime characterized by transfer of energy to high frequency 
modes [lj. We have uncovered that after this weakly nonlinear 
regime a turbulent regimes takes over. Crucially, this turbu- 
lent regime is characterized by a Kolmogorov-Zakharov scal- 
ing which we argue is compatible with 5/3. We have obtained 
similar results for the gravitational collapse in asymptotically 
AdSi spacetimes. 

We would like to highlight some of the key open problems 
that our results motivate. 

We have discussed that the steady-state assumption in tur- 
bulent flows is achieved in our problem through a separa- 
tion of scales. Given the connections between gravity and 
the Navier-Stokes equation [lj], it will be interesting to pur- 
sue a description of collapse more along the lines of turbu- 
lent fluid dynamics. Alternatively, we have indications that 
the phenomena we are witnessing in gravitational collapse is 
more along the lines of wave turbulence Ql| . The wave 
turbulence approach, being more general, might provide an 



explanation for various aspects that are not typical of full- 
blown Kolmogorov turbulence such as the similarity of our 
results in AdSs and AdSi. More importantly, wave turbulence 
might provide a semi-anlytical approach to our Kolmogorov- 
Zakharov spectrum. Aspects of wave turbulence in gravita- 
tional collapse will be reported elsewhere [T^. Similar appli- 
cations of wave turbulence have been developed for relativistic 
system [Isl.[l9j and provide a tentative explanation of why we 
find a Kolmogorov-Zakharov (rather than pure Kolmogorov) 
spectrum in a relativistic system. 

There is gathering evidence in the community of numerical 
and mathematical relativity that the product of the evolution 
studied in this paper which is a black hole is further unsta- 
ble. It would be very interesting to use similar methods to 
study the evolution of perturbations on top of a background 
that contains a black hole. Similarly, it has recently been 
noted that boson stars might be stable in asymptotically AdS 
spacetimes because the eigenmodes fail to meet the resonance 
condition [2(| • It would be interesting to scrutinize this claim 
under the non-linear perturbation microscope. 

Finally, given that the process of gravitational collapse is 
dual to thermalization in field theory, it makes sense to spec- 
ulate that the universality in the turbulent regime means that 
for strongly coupled field theories the process of thermaliza- 
tion always goes through a similarly turbulent channel. 

Acknowledgments We are thankful to C. A. Terrero-Escalante 
for collaboration on related topics and to D. Minic and D. 
Reichmann for important comments and insights. H. P. O. is 
grateful to the financial support of Brazilian agencies CNPq 
and PAPER J. E. L. R. thanks FAPERJ/CAPES for financial 
support. L.A.P.Z. is thankful to KITP and Aspen Center for 
Physics for hospitality during various stages of this work. This 
research was supported in part by the National Science Foun- 
dation under Grant No. NSF PHY11-25915 (KITP), grant 
No. 1066293 (Aspen) and by Department of Energy under 
grant DE-FG02-95ER40899 to the University of Michigan. 



[1] P. Bizon and A. Rostworowski, "On weakly turbulent 
instability of anti-de Sitter space," Phys. Rev. Lett. 107 
(2011) 031102 [arXiv: 1104.37021 [gr-acj], 

[2] J. Jalmuzna, A. Rostworowski and P. Bizon, "A comment 
on AdS collapse of a scalar field in higher dimensions," 
Phys. Rev. D 84 (2011) 085021 arXiv:1108.4539l [gr-qc]]. 

[3] O. J. C. Dias, G. T. Horowitz and J. E. Santos, "Grav- 
itational Turbulent Instability of Anti-de Sitter Space," 
larXiv:1109. 18251 [hep-th] . 

[4] D. Garfinkle, L. A. Pando Zayas and D. Reichmann, "On 
Field Theory Thermalization from Gravitational Col- 
lapse," JHEP 1202 (2012) 119 |arXiv:1110.5823l [hep- 
th]]. 

[5] D. Garfinkle and L. A. Pando Zayas, "Rapid Thermaliza- 
tion in Field Theory from Gravitational Collapse," Phys. 
Rev. D 84 (2011) 066006 [arXiv: 1106. 23391 [hep-th]]. 

[6] H. P. de Oliveira, L. A. Pando Zayas and C. A. Terrero- 
Escalante, "Turbulence and Chaos in Anti-de-Sitter 
Gravity," larXiv:1205.3232l [hep-th]. 

[7] A. N. Kolmogorov, "Local structure of turbulence in an 
incompressible fluid for very large Reynolds numbers," 
Doklady Acad Sci. USSR 31 (1941) 301-305. 

[8] S. B. Pope, "Turbulent Flows," Cambridge University 



5 



Press, 2000. 

[9] U. Frisch, "Turbulence: The Legacy of A. N. Kol- 
mogorov," Cambridge University Press, 1995. 

[10] J. Boyd, Chebyshev and Fourier spectral methods, 2nd 
edn. (Dover, New York, 2001). 

[11] A. Ishibashi and R. M. wald, Class. Quant. Grav. 21, 
2981 (2004). 

[12] S. Newhouse, D. Ruelle, and F. Takens, Comm. Math. 

Phys., 64, 35 (1978). 
[13] A. E. Motter, "Relativistic chaos is coordinate invariant," 

Phys. Rev. Lett. 91 (2003) 231101 [gr-qc/0305020] . 
[14] I. Bredberg, C. Keeler, V. Lysov and A. Strominger, 

"From Navier-St okes T o Einstein," JHEP 1207, 146 

(2012) [arXivillOl. 24511 [hep-th]]. 



[15] V. E. Zakharov, V. S. LVov and G. Falkovich, Kol- 

mogorov Spectra of Turbulence I - Wave Turbulence, 

Springer, 1992. 
[16] S. Nazarenko, "Wave Turbulence," Springer, 2011. 
[17] G. During, C. Falcon and L. A. Pando Zayas, " The Wave 

Turbulence Path to Gravitational Collapse," In progress. 
[18] R. Micha and I. I. Tkachev, "Relativistic turbulence: A 

Long way from preheating to equilibrium," Phys. Rev. 

Lett. 90 (2003) 121301 [hep-ph/0210202] . 
[19] R. Micha and I. I. Tkachev, "Turbulent thermalization," 

Phys. Rev. D 70 (2004) 043538 [hep-ph/040310l] . 
[20] O. J. C. Dias, G. T. Horowitz, D. Marolf and J. E. Santos, 

"On the Nonlinear Stability of Asymptotically Anti-de 

Sitter Solutions." TarXiv: 1208.57721 [gr-qc]. 



